REGULARITY AND DECAY OF SOLUTIONS OF NONLINEAR 

HARMONIC OSCILLATORS 



MARCO CAPPIELLO AND FABIO NICOLA 

Abstract. We prove sharp analytic regularity and decay at infinity of solu- 
tions of variable coefficients nonlinear harmonic oscillators. Namely, we show 
holomorphic extension to a sector in the complex domain, with a corresponding 
Gaussian decay, according to the basic properties of the Hermite functions in 
R''. Our results apply, in particular, to nonlinear eigenvalue problems for the 
harmonic oscillator associated to a real-analytic scattering, or asymptotically 
conic, metric in R'^, as well as to certain perturbations of the classical harmonic 
oscillator. 



1. Introduction 

The harmonic oscillator H = —A -|- in M'^ represents one of the simplest and 
yet more useful models for several physical phenomena, and its relevance both in 
Mathematical Analysis and Physics is well-known. Its eigenf unctions, namely the 
Hermite functions ha{x), are given by the formulae ha{x) = pQ,(x)e~'^'' a G N"', 
where pa is a polynomial of degree \a\ (see e.g. [2H])- Two remarkable features 
of the Hermite functions are their Gaussian decay at infinity, and their very high 
regularity. In fact, we have 

(1.1) |/ia(a;)| <e-'=l^l' for x G K(OI ^ e"^'^'' ioi i eW^ 

for every c < 1/2, where ha{C) denotes the Fourier transform of ha- The functions 
ha in fact extend to entire functions ha{x + iy) in the complex space and, for 
every < e < 1, we have the estimates 

(1.2) \ha{x + iy)\ < e'"^""^^ in the sector |y| < £:(1 |a;|), 
for some c > 0. 

In this paper we wonder to what extent these properties continue to hold for non- 
linear perturbations of the harmonic oscillator, possibly with variable coefficients. 
Relevant models are equations of the type 

(1.3) - Au + \x\\ - Xu = F[u], A G C, 
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with a nonlinearity of the form F[u\ = X]|q:|+|/3|<i Ca/sa^^f^'^w'', k > 2. Cappiello, 
Gramchev and Rodino in [7J showed by a counterexample that generally, even in 
dimension d = 1, there can exist Schwartz solutions of f ll.Sp which do not extend 
to entire functions in C. In fact, a refinement of their argument (see Section [5] 
below) shows that a sequence of complex singularities may occur, approaching a 
straight line at infinity. On the other hand, as a positive result, it was proved in 
[7] that every solution u G H^{M.'^), s > d/2 + 1, of (\1.3\i extends to a holomorphic 
function u{x + iy) on the strip {2; G : \lmz\ < T} and satisfies there an estimate 
of the type \u{x + iy)\ < Ce~^^^^ , for some c,C,T > 0. Similar results, namely, 
holomorphic extension to a strip and super-exponential decay, were proved in [71 [10] 
for more general classes of elliptic operators with polynomial coefficients. 

The above mentioned negative result as well as the estimates f ll.2p . valid in a 
sector in the linear case, suggest the possibility, even in the presence of certain non- 
linear perturbations, of a holomorphic extension of the solutions to a sector, rather 
than only a strip, with a corresponding Gaussian decay estimate. In this paper 
we show, for a large class of equations including (11. 3p . even with non-polynomial 
coefficients, that this is in fact the case. The techniques developed here actually 
will apply to much more general differential (and pseudodifferential) operators. 
To motivate the class of operators we will consider, we first discuss a special yet 
important example. 

Consider the equation Pu = F[u], with 

d d 

(1.4) P=Y1 9'\x)d,dk + ^k{x)dk + V{x), 

j,k=l k=l 

where the functions g^'^, bk, and the potential V are real- analytic in R*^, and satisfy 
the following conditions. 

We suppose that the matrix [g-''') is real and symmetric and that there exists a 
constant C > such that 

d 

(1.5) 5Z^7^''(x)e,e.>c7-i|eP vx,eeM^ 

j,k=l 

as well as 

(1.6) \d^g^''{x)\ + \d^bk{x)\ < Va; G a G 
where (x) = (1 + |xp)"^/^. Moreover we assume that 

{ReV{x)>C~^\x\'^ for|x|>C, 
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We consider a nonlinearity of the form 

(1.8) F[u]= Fju{x)u\dju)\ 

2<h+l<N 
\<i<d 

for some G N, where 

(1.9) \d''FjM{x)\ < Cl°l+ia!(a;)i""^^^^i''>-l"l Vx G a G N'^. 

Then, we claim that 

Under these assumptions, every solution u G if^(M°'),s > d/2 + min{l,/}, of 
the equation Pu = F[u], extends to a holomorphic function u{x + iy) in the sector 
{z = X + iy E : \y\ < s{l + \x\)} of for some e > 0, satisfying there the 
estimates \u{x + iy)\ < Ce~^'^' , for some constants C > 0, c > 0. 

Notice that if V{x) satisfies (11.71) then also V{x) — A, A G C, satisfies it, so that 
the above result applies to the corresponding eigenvalue problem as well. In the 
linear case {F[u] = 0) this result intersects the wide literature on the decay and 
regularity of eigenfunctions of Schrodinger operators, cf. Agmon [1], Nakamura 
[22] . Sordoni [27], Rabinovich [21], Rabinovich and Roch [22] and many others. 

We also remark that suitable perturbations of the standard harmonic oscillator 
fall in this class of equations, as well as the harmonic oscillator associated to a 
real-analytic scattering, or asymptotically conic, Riemannian metric in M.'^ (see 
Section [5] below) . For a detailed analysis of these metrics and their important role 
in geometric scattering theory we refer to Melrose [T9l [20] , Melrose and Zworski 

m- 

Let us now state our main result in full generality. We consider nonlinear equations 
whose linear part is a differential or even pseudodifferential operator 

(1.10) Pu{x)=p{x,D)u{x) = {2'k)-^ j e'''^p{x,i)u{i)di, 

with symbol p in the class T^(W^),m > 0, defined as the space of all functions 
p G C°^(]R^'^) satisfying the estimates 

(1.11) \dp^p{x,0\ < Cl"l+l^l+^a!/3!(l + |a;| + |e|)'"-l"l(a;)-l^l 

for all (x, e M^^, a,/3 e N"', and for some positive constant C independent of 
a, (3. This class is particularly suited to study harmonic oscillators with variable 
analytic coefficients and it is inspired by the class considered by Shubin [26] and 
Helffer [15] which was in fact modelled on the harmonic oscillator and its real 
powers. However, a differential operator belongs to that class only if its coefficients 
are polynomial, which is an unpleasant limitation. With respect to [15], [26], we 
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overcome this restriction by assuming the less demanding estimates (11. lip , which 
at the same time imply that the symbol p is real-analytic. For example, a function 

\a\<m 

belongs to r^(M'^) if the coefficients c„ satisfy |<9''c„(a;)| < Cl''l+i/3!(x)'"-l"l-l'^l. 
We shall assume moreover the symbol p of our operator to be F-elliptic, in the 
sense that, for some constant R> 0, 

(1.12) inf (l + |x| + |e|)-™b(x,OI>0. 

\x\ + \^\>R 

This is clearly a global version of the classical notion of ellipticity. For example, 
by (ll.5p -( lL71) . the symbol of the operator in (II. 4p belongs to F^(M'^) and satisfies 
( irT2|) with m = 2. 
We moreover consider a nonlinearity of the form 

(1-13) F[u]= Yl F,U.-MU^'"^, 

h,l,pi,...,pi k=l 

where the above sum is finite and h,l ^ N,l > 2, pi, . . . , pi & N'^ satisfy the condi- 
tion h + max{|pfc|} < max{m — 1,0}. We assume that the functions F^ip^ p^^x) 
satisfy the following estimates 

(1.14) |9''F,,,,,,...,,(x)| <Cl^l+i/3!(x)'^-l^l, 

for some positive constant C depending on h,l, pi, . . . , pi and independent of (3. 
Our main result is the following. 

Theorem 1.1. Let p G F™(M'^),m > 0, satisfy (I1.12p and let F[u] be of the form 
(ll.lSp . (I1.14P (possibly with some factors in the product replaced by their conju- 
gates). Assume that u G H^{W^), s > (i/2 + maxfc{|pfc|} is a solution of the equation 
Pu = F[u]. Then u extends to a holomorphic function u{x + iy) in the sector 

(1.15) {z = x + iy eC^ : \y\ < e{l + \x\)} 
of C^, for some e > 0, satisfying there the estimates 

(1.16) \u{x + iy)\<Ce~'^'^^\ 

for some constants C > 0, c> 0. The same holds for all the derivatives of u. 

The linear case F[u] =0 deserves a special interest and will be treated in detail 
in Section |5l We emphasize the fact that the form of the domain of holomorphic 
extension as a sector is, in a sense, completely sharp, even for the model (II. 3p (see 
Section [5]) . 
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Let us briefly compare our result with those in the existent literature. Several 
papers were devoted to the problem of holomorphic extension to a strip and ex- 
ponential decay of solutions of certain semilinear elliptic equations arising in the 
theory of solitary waves or bound states, whose model is —Am + u = \u\p~^u, 
cf. Beresticky and Lions [2j, Bona and Grujic' [1], Bona and Li [5], [6], Biagioni 
and Gramchev [3j, Gramchev [13], Cappiello, Gramchev and Rodino [8], [9] and 
the references therein; see also our recent paper [TT] for the extension to a sec- 
tor. However, as it is clear from our model fll.3p . we consider a different class of 
equations here, and in fact we deal with Gaussian, rather than exponential decay. 
Instead, as already mentioned, a class similar to the present one was considered in 
[7], where the problem of the extension to a strip, combined with super-exponential 
decay, was addressed. The main novelties of the present work are the possibility 
of treating non-polynomial coefficients and nonlinearities, and the achievement of 
the optimal extension result, namely to a sector. 

The paper is organized as follows. In section [2] we list some known factorial and 
binomial estimates and we collect some basic properties of the pseudodifferential 
operators introduced before, which will be instrumental in the proofs of our results. 
In Section E] we introduce a suitable space of analytic functions which exploits the 
two properties fll.lSp and fll.l6p . Section H] is devoted to the proof of Theorem ll.il 
which is based on an iterative scheme on the space defined in Section [31 Finally, in 
Section [5] we give some concluding remarks. In particular, we read our results on 
the models introduced above and treat in detail their application to the Schrodinger 
operator in M*^ with a scattering metric. Finally, we discuss the sharpness of our 
results for what concerns the shape of the domain of the holomorphic extension as 
a sector of C^. 



2.1. Factorial and binomial coefficients. We use the usual multi-index nota- 
tion for factorial and binomial coefficients. Hence, for a = (ai, . . . , a^) G N"^ we 
set a! = ai! . . . a^! and for (3,a E N'^, (3 < a, we set (j^) = °_|_^^| . 
The following inequality is standard and used often in the sequel: 



2. Notation and preliminary results 



(2.1) 




Also, we recall the identity 




j = 0, l,...,|a|. 
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which follows from ULii^ + ^T" = (1 + ^)'°'' and gives in particular 

(2.2) G)^Q' ".^^«^/'^- 

The last estimate implies in turn, by induction, 

a\ \a\\ 
5i\...5j\ - \5i\\ . . .\5j\\ 

as well as 



(2-3) y^-TT< I.I, I. IP « = 5i + ... + 5„ 



Finally we recall the so-called inverse Leibniz' formula: 

(2.5) xPd''u{x)= {^^d"~\^^~^u{x)). 

2.2. Pseudodifferential Operators. We collect here some basic properties of the 
class r™(R°') defined by the estimates (11.111) and of the corresponding operators 
(ll.lOp . Actually, for our purposes it is not necessary to develop a specific calculus 
for the analytic symbols of r™(]R'^). We shall deduce the properties we need from 
those of the larger class r™(]R'^) defined as the space of all functions p e C°°{M?'^) 
satisfying the following estimates: for every G N'^ there exists a constant 
Ca,i3 > such that 

(2.6) \d!d^p{x,0\<C.A^ + \x\ + \C\r-\'^Hx)-\^^ 

for every x,^ G R'^. Clearly r™(M'^) C r"(M'^). We shall denote by OPr"^(R'^) 
(respectively OPr^(M'^)) the class of pseudodifferential operators with symbol in 
r"^(M'^) (respectively in r^(M'^)). We endow r™(M'^) with the topology defined by 
the seminorms 

\\ptP = sup sup {|95"9fp(x,OI(l + kl + |el)"'"+'"'(a;)l''l}, iVGN. 

\a\ + \l3\<N (a:,^)eM2d 

The properties of r'"(M'^) follow from the general Weyl-Hormander's calculus in 
[T71 Chapter XVIII]; with the notation used there, r™(M'^) = S{M,g), with the 
weight M(x, ^) = (1 + |xp + I^P)""/^ and the metric 

. ^ \dy\' , \dv\' 

^-^^^'^^ = rTN^ + i + ia:p + ier 

We also refer the reader to |23l Chapter 1] for an elementary and self-contained pre- 
sentation; with the notation in [231 Definition 1.1.1] we have r'"(M'^) = S{M] $, '^), 
with M{x,C,) as above and 

$(x,o = (x), ^{x,o = a + \x\' + \^\Y'- 
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Now, if p G r'"(]R'^) then p(x,D) defines a continuous map 5(R'^) — )• 5(M'^) 
whicli extends to a continuous map — )■ iS'(]R'^). The composition of two such 
operators is therefore well defined in 5(M'^) and in iS'(]R'^); more precisely, if pi G 
r"^i(M'^) andp2 e T"'^(R'^),then pi{x, D)p2lx, D) = p^ix, D) with p^ G r'^i+'^2(M^) 
and the map (^1,^2) ^ P3 is continuous r™i(M'^) x r™2(M<i) r™3(R'^). Moreover 
we have that f] r™(]R'^) = iS(R^'*). In particular, operators with Schwartz symbols 

are globally regularizing, i.e. they map continuously 5'(M'^) into 5(]R'^). 

Operators in OPr™(]R'^) are also bounded on certain weighted Sobolev spaces. 
We consider, for simplicity, the case of integer positive exponents (we will only 
need this case, see [2S] for the general case). For s G N, we define 

(2.7) Q"(M'^) = {ue L^R"^) : ||m||q. := ^ \\x^d''u\\L2 < 00}. 

We recall that f] Q"(M^)>S(M^). Now, if P G OPr'^(M'^), m G Z, m < s, we have 

sgN 

P : Q'(R'^) Q"-"^(M^) continuously with 

|b(x,D)|U(Q.,Q.-™)<Clb||P 

for suitable C > 0, G N depending only on s, m and on the dimension d; see [231 
Proposition 1.5. 5, Theorem 2.1.12]. Moreover, for m & Z, m < s, there exists an 
operator T G OPr-™(]R<^) which gives an isomorphism Q"-™(]R'^) Q'iR'^). 

We will also need the following Schauder's estimates for the weighted Sobolev 
spaces Q'(R'^) in (lO) . 

Proposition 2.1. Let s G N, s > d/2. There exists Cg > such that 

\\uv\\qs < Cs\\u\\qs\\v\\qs \/u,V G (5^(M'^). 

Proof. We have 



\uv\\qs= J2 \\x^d''{uv)\\L2= ^l^yix^d^-^u-d^vW 



L2 



q| + |/3|<s |a| + |/3|<s 7<" 

q:| + |/3|<s 7<a 



<2^ V ^Wx^d^-^uWLvWvWL-, 



where I < p, q < 00 are chosen to satisfy - + - = h, and i<- + 4^, i<- + 



p ' q 2' 2 p ' d ^ 2 q ' d 



This is possibile because s > f . Then, by the Sobolev embeddings we have 



(2.8) \\uv\\qs<Cs Y1 J2\\^'^^"~^'^\\Hh\\\d^v\\H.-i-,\ 

\a\ + \l3\<s 7<a 
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On the other hand, 

(2.9) X J2 ll^'' (x^d^-^u) ||i2 < 

ImI<|7| 

Similarly, 

(2.10) WvWh^^n < C':\\u\\qs. 

Combining fl2.8p . fl2.9p and f l2.10p we get the desired result. □ 
As a technical tool, we will also use the scale of weighted Sobolev spaces 

(2.11) H'^'^^R'^) = {ue 5'(M^) : Mh^s^ := || (a;)"^^!!,, < oo}, 

defined for Si,S2 G M. In particular, we need the following result (see e.g. [221 
Definition 3.1.1 and Theorem 3.1.5]). 

Proposition 2.2. Consider a symbol p{x,^) satisfying the estimates 

(2.12) \d^d^p{x,0\ < Ca,/3(a;)""'^'(0""'"', Va,/3 G N'', x G R''. 

Then the corresponding operator p{x, D) is bounded H'^^''^'^{R'^) — )■ H''^'^~"^'''^^''^(R'^) 
for every Si,S2 € R, with operator norm estimated by an upper bound of a finite 
number of the constants Ca^j3 appearing in fl2.12p . 

By using Schauder's estimates in the standard Sobolev spaces and the inclusion 
H^i^'^iW^) ^ ff^i(M^), valid if S2 > 0, one also gets 

(2.13) ||ut;||Hn,s2 < Csi,s2ll^||//"i'=2 ||t'||ii'^ii.«2 si > -, S2> 0. 

A symbol p G r'^(R'^) (and the corresponding operator) is said to be T-elliptic if 
it satisfies the condition 01.120 . 

The notion of F-ellipticity for an operator in OPr™(R"') will be crucial in the 
subsequent arguments because it guaranties the existence of a parametrix E G 
OPr~'"(R'^). Namely we have the following result, see [23| Theorem 1.3.6] for the 
proof. 

Proposition 2.3. Let p G r'"(]R'^) be T-elliptic. Then there exists an operator 
E G OPr-™(R'^) such that EP = I + R and PE = I + R' , where R, R' are 
globally regularizing pseudodifferential operators, i.e. R and R' are continuous 
maps S'{W^) — ^ iS(R''). The operator E is said to be a parametrix for P. 

Finally we point out for further reference the following formulae, which can be 
verified by a direct computation: for a, /3 G N*^, u G iS(M'^) : 

(2.14) x^p{x, D)u = V(-l)l^l {Djp){x, D){x^-'^u), 



7</3 ^ 
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(2.15) crp{x,D)u = Y,(fj 



idip){x,D)d"~^u. 



3. A SPACE OF ANALYTIC FUNCTIONS 

We introduce a space of real-analytic functions in M*^, which extend holomorphi- 
cally on a sector in and display there a Gaussian decay. 

Definition 3.1. We denote by nsect(^'^) the space of all functions f E C°°(]R'^) 
satisfying the following condition: there exists a constant C > such that 

(3.1) |x^9"/(x)| < Cl"l+l^l+^M(a,/3), for all a, /3 G 
where 

(3.2) M{a, (3) = inax{|a|, 

It is easy to verify that the space 'HsectiM'^) is closed under differentiation. 

Theorem 3.2. Let f G 'Hsecti^'^)- Then f extends to a holomorphic function 
f{x + iy) in the sector 

(3.3) Ce = {z = x + iy eC^ : \y\ < e{l + \x\)} 
of for some e > 0, satisfying there the estimates 

(3.4) \f{x + zy)\<Ce-^\^\\ 
for some constants C > 0, c > 0. 

Proof. First we show the estimates 

(3.5) |x^9"/(x)| < Cl"l+>|!e-'^l^l', for \(3\ < \a\. 

Indeed, since < fc"X]|7|=n l^^''^! ^ constant k > depending only on the 
dimension d, by (13. ip we have (assuming C > 1 in (13. ip ) 

^ (clxP''" 



e-\-\'\x^d"f{x)\ = J2 , l^^'9"/(x)| 



ni 

n=0 

oo ^ 

<Y.(^krY.—y^'wy{x)\ 

n=0 |7|=n 
oo 

n=0 h\=n 

CO 

n=0 |7|='^ 





|!i/2(|« 


| + 2|7|)!^/^ 




7l 


I 
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where in the last step we used the inequahty {\a\ + 2|7|)! < 2l"l+^l'''l|a|!|7|!^, which 
follows by applying twice f l2.ip . Since the number of multi-indices 7 satisfying 
I7I = n does not exceed 2'^+""^, we get (13 .Sp for a new constant C, if c is small 
enough. Now, (13. 5p and the estimate |a|! < c/l^'a! give 

(3.6) |5"/(x)| < Cl"l+i«!(x)'l"le-"l^l', 

for a new constant C > 0. By considering the Taylor expansion of the function 
/ centered in any x G M*^ and using the estimates in (13. 6 p we obtain the desired 
extension property in a sector of the type (13.31) together with the estimates (13. 4p . 
□ 

In the sequel we will use the following characterization of the space HsectiJ^'^) in 
terms of Q'^-based norms. 
Set, for / e 5'(R"'), 

where M(a,/3) is defined in (13. 2p . 

Proposition 3.3. Let f G 'Hsecti^'^)- Then for every s G N there exists e > such 
that S:;,'[f]< 00. 

In the opposite direction, if for some s G N there exists e > such that S^[f] < 
00, then / G "H. 



ecty 



Proof. Assume / G 'Hgecty^ )■ We have 

|5| + |7|<s 

Now, if M G N satisfies M > d/Awe have 

(3.8) \\x^d^{x'^d''f)\\L2 < C"||(l + |a;n^V9^(x^9"/)|Uoo. 

By Leibniz' formula, (13. ip and (12. ip we get 

||x^5"/||q=<cI"|+i^i+'m(«,/3) 

for some constant > 0. Hence S^[f] < 00 if £ < C^^. 

In the opposite direction, we may take s = 0; hence assume < 00 for 

some e>0. Then ||a;^(9"/(a;)||L2 < Cl"l+l^l+iM(a, /3) for all a, (3 e N"'. If M is an 
integer, M > d/2, we have 

\\x^d'-f\\L^<C J2 W{x^d"f)\\L2, 

h\<M 

and similarly one gets that / G Hsecti^'^)- □ 
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4. Proof of the main result (Theorem 11.11 ) 

In this section we prove Theorem 11.11 In fact we shall state and prove this result 
for the more general non-homogeneous equation 

(4.1) Pu = f + F[u], 

where P and F[u] satisfy the assumptions of Theorem II. II and / is a function in the 
space 'Hsecti^'^) defined in Section 3. Moreover we can restate our result in terms 
of estimates in 'Hsecti^'^)- Namely, in view of Theorem 13. 2 1 it will be sufficient to 
prove the following theorem. 

Theorem 4.1. Let P = p{x, D) G OPrj^(M'^), m > 0, he T-elliptic, that is ffLT2|) is 
satisfied. Let F[u] be of the form f ll.l3p (possibly with some factors in the product 
replaced by their conjugates) and f G 'Hsecti^'^)- Assume moreover that u G H'^iW''), 
s > d/2 + maxfc{|pfc|}, is a solution of (14. ip . Then u G 'Hsect{W^)- 

In fact we always assume that F[u\ has the form in (11.131) . and we leave to the 
reader the easy changes when some factors of the product in (I1.13P are replaced by 
their conjugates. 

The first step is to show that, under the assumptions of Theorem 14. H the sum 
S'^^[m] is finite for every G N and for some e > 0. In particular, we prove the 
following preliminary result. 



Lemma 4.2. Under the assumptions of Theorem 4-1, we have u E S{ 



Proof. The proof is based on a bootstrap argument in the scale of weighted Sobolev 

spaces defined in fl2ll]) . Notice that H'(R'^) = H'^^(R'^) and that f] H'^^'^iR'^) 

si,s2eR 

= iS(M'^). To prove the lemma it is then sufficient to show that if w G if^i'*2(M'^), 
si > d/2 + M, with M = maxk{\pk\}, S2 > 0, then u G i7"i+^.«2+^, with r = 
min{m/2, 1/2}, and to iterate this argument. We first consider the case m > 1. 
Let E G OPr"'"(R"') be a parametrix of P (Proposition EJ]). Applying E to both 
sides of the equation (14. ip , we get 

u = -Ru + Ef + EF[u], 

where R is globally regularizing. In particular, we have Ru G 5(M'^) and Ef G 
5(R'^) because / G 5(]R'^). Concerning the nonlinear term, we observe that, since 
—m < —M — h — 1, then the symbol e{x,C,) of E satisfies the following estimates 

\d^d^Mx,0\ < C,/3(l + |x| + |e|)^'^-'-^^"l(a:)-l^l 
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It follows from Propositfon O that E : H'^~^^,^2-h^jgd^ ^ jjs^+i/2,s2+i/2^-^d^ 
tinuously for every Si, S2 G M. In particular, in view of (11 .Mp . we have 



+ 1/2,32 + 1/2 



\\EF[u]\\hsi+i/2,s2+i/2 = \\E Fh,i,p,...pi{x)Y[d'"'u\\Hsi 

h,l,pi,...,pi k=l 

I 

h,l,pi,...,pi k=l 

I 

< C'\\Y[dP'u\\Hs,-M.s, < C"\\ 

k=l 

by Schauder's estimates fl2.13p . because Si — M > d/2. The case < m < 1 is 
completely similar, considering that in this case h = M = and the symbol e{x,C,) 
of E satisfies the estimates 

so that E maps continuously H''''^{R'^) into /7«i+™/2,s2+m/2(^K'^). □ 

In order to prove Theorem 14.11 it suffices to verify that [u] < 00 for some 
s > 0,£ > 0, in view of Proposition 13.31 This will be achieved by an iteration 
argument involving the partial sums of the series in (13.71) . that are 

(4.2) s-'if]= ^ ^!^||.fa«/i|,., 

\a\ + \l3\<N '^-'V"./^; 

where M(a,/3) is defined in ( 13. 2p . 

4.1. Proof of Theorem 14. IL We need several estimates to which we address now. 

Proposition 4.3. Let R G OPr^^(]R'^). Then for every s G M there exists a 
constant Cg > such that, for every e > 0, N > 1 and u G 5(]R'^), we have 

o<H+W<N ^ 

Proof. We first estimate the terms with /3 = 0, hence a 7^ 0. Let k G {1, . . . , c?} 
such that afc 7^ 0. Since RodkE OPr°(M'^) is bounded on Q"(M^) we hav(l3 



We denote by Ck the fcth vector of the standard basis of M . 
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On the other hand, when /3 7^ 0, hence f3j 7^ for some j G {1, . . . ,d}, we use 

j 



the fact that Roxj E OPr°(R'^) is bounded on Q"(M^). We get 



Since M(a, /3) > M(a, {3 — ej), this gives the desired result. □ 



Proposition 4.4. Let P = p{x,D) be a pseudodifferential with symbol p{x,C,) 
satisfying the estimates (11. Ill) , with m > 0. Let E G OPr~™'(]R'^). Then for every 
s G M there exists a constant Cs > such that, for every e small enough, N > 1 
and u G 5(]R'^), we have 

(4.3) M7^"^f^''^'^"^''"«^ - ^^"^^-iM- 

Proof. We estimate separately each term arising in the sum (14. 3p . We write 

[p, x^d""] = [P, x'^ja" + x^[P, d"]. 
Hence, using (I2.14p . (I2.15p . we get 

(4.4) [P,x^d'']u= y i-iy-'^'^+^f ^\D]'p){x,D){x^-^'d''u) 

(^\xf'dip{x,D)cr-^u. 



07^70 </3 



Given /3, 5, let 5 be a multi-index of maximal length among those satisfying \5\ < \5\ 
and 5 < (5 (hence, if \5\ < \5\ then /3 — 5 = 0). Writing = x^x^~^ in the last 
term of (14. 4p and using again (I2.14p we get 
(4.5) 

[P,x'^9"]m = Y Y1 (-1)'^°'+'(^ ~ ^) h\x\Dj''d^,p){x,D){x^-~^-^'>d''-^u). 

(5,70)^(0,0) 

We now work out this formula to obtain a useful representation of the commu- 
tator [P, Namely, we look at the operator x^~^~'^°d"^~^ . Given 70, a, 5, let 
7o be a multi-index, to be chosen later on, satisfying |7o| < |7o| and 70 < a — 5. 
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We write, by the inverse Leibniz formula fl2.5p . 

0^ (/3-^~-7o-7i)! V7i; 

0/71</5-i5 — I'D 
7l<70 

We now look at the operator x^~^~^"~'^'^d°'~^~^° . We denote by 71 a multi- index, 
to be chosen later on, satisfying |7i| < |7i|, 71 < a — 6 — %. Applying again the 
inverse Leibniz formula we have 

(4 7) d"^^^^'^ = x^~^~^°~'^'^ d^^ d°'~^~^°~^^ 

_j_ \ ^ ( — 1)1^ ^(/^ ~ 1^ ~ 70 ~ 7l)- f 7l\ q')i-12 q ^;S-5-70-7i-72 ^a-<5-7o-7i 

0^ Si 5-70 -71 -72)! V72/ 

72<71 

Continuing in this way and substituting all in (14 .Sp we get 

r 

(4.8) [P,a:^a"]« = EE E E E ^"./3A7o,7....,7. 

<5<« i=0 7o<,9-6 0#7i</9-«-7o 0#7j </3-i-7o-...-7j_l 
(5,70)^(0,0) 7i<70 7j<7j-i 

XPa,/3,5,7o,7„...,7. (a;, ^) , 

where 7^ satisfy |7j| < |7j| and 7j < a — 5 — 70 — ... — 7j-i, 

(4.9) P„,;3,5,70,7l,...,7.(^,O = J > 0, 

and 

' (« - 5)!5!7o!(/5 - 5 - 70 - ... - 7.)! M V 7. 

(4.10) < l"^!'!^^' '^1- 2lT'o+-+'^^-il, 

|a - 5|!5!7o!|/3 - 5 - 7o - . . . - 7j|! 

cf. (12. 4p and (12. ip . (If j = in (I4.10p we mean that there are not the binomial 
factors, nor the power of 2). Observe that, since we have 7^ 7^ for every j > 1, 
this procedure in fact stops after a finite number r of steps. 

We now study separately the terms with \a\ > \(3\ and \/3\ > \a\. 

Case \a\ > 1/31. 
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We use the formula f l4.8p . where now we choose 70 satisfying, in addition, |7o| = 
|7o|. Such a multi- index exists, because |a| > \/3\. Similarly, at each subsequent 
step we can choose 7j, j > 1, satisfying in addition = |7j|. 

Now we observe that, by (11. lip . (12. ip . and Leibniz' formula, for every 6,a E N'^ 
we have 

(4.11) |9|9:p„,^,5,,o,7i,...,7.(^,oi < ci^°i+i^i+So!<5!(i + \x\ + 

for some constant C depending only on 9 and a. In fact, \6\ < \6\, \% — 71 + 
7i — . . . — 7j + 7j| = |7o| = I70I, and the powers of \S\ and |7o| which arise can be 
estimated by C^'yo\+\s\+-^ for ^.q^q C > 0. 

We now use these last bounds to estimate E o Pa,p,s,'yo,ii,-nA^^ ■^)- ^o this 
end, observe that this operator belongs to OPr°(M'^), and therefore its norm as a 
bounded operator on Q'*(M'^) is estimated by a seminorm of its symbol in r°(M'^), 
depending only on s and d. Such a seminorm is in turn estimated by the product 
of a seminorm of the symbol of E in r~™'(M'^) and a seminorm of Pa,/3,(5,7o,7i,...,7j in 
r™(M'^), again depending only on s,d. Hence, from (14. lip we get 

(4.12) \\Eop^^p^,^^^^^^_^^{x,D)\\,sm < C^^I+l'l+So!5!. 
Since |7fc| = |7a;|, < A; < j, we have 

(4.13) \P-ma-5-^,-...-^,\\ ^ ^ 

\a - 6\\\I3 - 5 - 7o - . . . - 7j|! 

(recall that if \5\ < \5\ then /3 — 5 = 70 = . . . = 7^ = 70 = . . . = 7j = 0). 
By fOUp . (021) Km . we get in this case 



(4.14) 

1/31 

\Co,,p,5,^,,^,_^^\\\Eop^,^^Sa^^^^^^ 

^|a| + |/3|-|<5|-|<5|-|7o+...+7j|-|7o + ...+7i 



<5| + |<5| + |7o+---+73l + l7o + ---+7j 



|a - 5 - 7o - . . . - 7j|! 

X 1 1 ^/3-5-70 - . . . -7. ^a-'5-70 - . . . -7. ^ 1 1 



Now, the assumption \a\ > and the choice of 5 and 7^, j > 0, imply M(a, (3) 
\a\\ and |/3 — 5 — 70 — . . . — 7j| < \a — 6 — % — . . . — Hence 

M{a - 5 - 7o - . . . - 7j, /3 - (5 - 7o - . . . - 7i) = |a - - 7o - • • • - 7i|!- 



u 
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We can therefore rewrite (14.141) as 
(4.15) 

< C'^(c^£:)I'5I+I^I+Ito+---+'^j1+It'o+---+^jI X 

^|Q| + |/3|-|<5|-|<5|-|7o+...+7j|-|7o+---+7jl 

M{a - 5 - 7o - . . . - 7j, /3 - 5 - 7o - . . . - 7j) 

We now perform the change of variables a = a — 6 — % — — jj, 13 = (3 — 5 — 
7o - . . . - 7j. In fact, the map (a, /3, 5, j, 70, 71, • • • , 7j) (a, /3, 5, J, 7o, 7i, • • • , 7j) 
defined in this way is not injective, because of the presence of 5, 70, . . . , 7j (of 
course, one should think of 5 as a function of a, /3, 5, and to every 7^, j > 0, as a 
function of a,f3,6,'yk, k < j). Anyhow, since \6\ < \6\ and = the number 
of pre-images of a given point is at most 2^^^~^^"'°^^'"^^"'^^~^'^^^~^'^\ Hence we deduce 
from (14. 8 p and (I4.15P that, if e is small enough, 

|a| + |;3|<JV ^ ^ ' 



|ci| + |^|<7V-l ^) i = 5 717^0,...,7,#0 

7o:(<5,7o)7^(0,0) 
r 

j=Q 

Case 1/3 1 > Here it is convenient to get separate estimates when |x| is large 
or small at the scale To make this precise, consider a function ip G C^(M'^), 

if{x) = 1 for |x| < 1 and if{x) = for \x\ > 2. Let then 

ipp{x) ^ 



(notice that /3 7^ 0, because of our hypothesis \f3\ > \a\). Hence v^^(x) = 1 for 



\x 



< and ip/six) = for \x\ > 2|/3|^/^. Moreover, we have 

(4.17) \d''vp{x)\ < 7 e X G 



for constants > 0. 



We write 

[P,x^d''] = ip^{x)[P,x'^d''] + (1 - ipp{x))[P,x^d° 
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and we split consequently the terms in f l4.3p . 

Estimate of ^^ \\E{{1 - ipp{x))\P,x^d'']u)\\Qs. 

We use the expression in (14. 5 p for [P, x'^S"], and we split the second sum, by 
considering separately the terms with — 5 — 7o| < |a; — 5| or — 5 — 7o| > |a; — 5|. 
This is equivalent to saying — 7o| < |«| or — 7o| > |a| because > \a\ implies 
\5\ = \6\. Moreover we apply the iterative argument at the beginning of the present 
proof to the terms with 1/3 — 70 1 < We obtain 

(4.18) M(^"^^^^ " MmP,^^d-H\Q^ < in + W 



where 



(<5,7o)^(0,0), |/3-7o|>|a| 

X \\E{x\l - ^p{x)){Dj'd^^p){x,D){x^-~^-^''d''-^u))\\Qs 



whereas 
(4.20) 

w-EE E E - E 

S<a J=0 70</3-*: 07^71 </3-i5-7o 07^7j </3-i5-70 - ■ • • ~7j _ 1 

(<5,7o)^(0,0), |/3-7ol<|a| 7i<70 7j<7i-i 

where the multi-indices 7^ will be chosen later on, satisfying < |7j| and 7^ < 
a — (5 — 7o — ... — 7j-i; the constants Ca,i3,5,-yo,'yi,...,^j satisfy (I4.10p . and (14.90 holds. 

Estimate of the terms in (I) (hence |/3 — 7o| > |«|) 

Since \6\ = \6\, by Leibniz' formula, (II. lip and (12. ip . for every 6,a eN'^ we have 

(4.21) \dld^{x\Dj°dip){x,0)\ < Cl^ol+l^l+So!'5!(l + kl + |e|)"^-l^°l-l^l(x)~l"l 

for some constant C depending only on 6 and a. 

Since 1 — y9^(x) is supported where > using Leibniz' formula again and 

flirrp we get 



dld:ix\l-Mn){Dj'^dtp)ix,0) 



< 



^|70| + |5| + 1^^!^!| 



l70l 



;i + i^i + iei)' 



(x) 
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for some constant C depending only on and a. As a consequence, 

(4.22) \\Eo{x\l - v,{x)){Dfdip){x,D))h^Qs^ < C^'+I^l+Sol^!!/?!"^. 
On the other hand, we have 

(4.23) {'-'M"^<^ 



7o 







I 


|/3-5-7o|!| 


a 







as well as 



(4.24) 



W'\aY^'^\P-~5-^,\\\a-5\\ 



^ — 

I <i 



\\a-5\m\J V|/3_5-^o|! 
^/ 

By diss), (H:^ and fOI|) we obtain 



„ , , ^|a| + |,3|-2|<5h|7o| 

< CJC,£)2|^l+l^o| ^ ^ ||x^~^-^°5"-^m|Iqs. 

" ^ ^ |«-(5|!V2|/3_5_^o|!i/2" '"^ 

Since > |q;| and |/3 - 7o| > then we have M(a,/3) = and 
M{a -S,f3- 6 -jo)\a - 6\\^/'^\f3 - 5 - 7o|!^/^ so that (K^ can be rephrased as 



^|a| + |/3|-2|6|-|7o| 



- ^ ^ M(«-5,/3-(5-7o) 
Estimate of the terms in (II) (hence \(3 — 7o| < |«|). 

In the iterative argument which led to (14.201) . we choose the multi- indices jj, j > 0, 
in the following way: 70 is a multi-index satisfying, in addition, |/3 — 5 — 7o| = 
\a — 6 — 7o|. Such a multi-index exists, because \6\ = \6\, \/3\ > \a\ and \f3 — 
7o| < moreover |7o| — |7o| = — \a\. Similarly, we can choose 71 satisfying 
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1/3 — 5 — 7o — 7i| = la — 5 — 7o — 7i|; in particular |7i| = |7i|. In general we can 
choose 7j such that 

(4.27) 1/3 - 5 - 70 - ... - 7,1 = |« - 5 - 70 - ••• - 7,1; 

hence = |7j| if j > 1. 

Notice that now in (14.91) we have \5\ = \6\ and |7o — 7i + 7i — • • • —7, +7jl = l7o|- 
Hence, since |7o| — |7o| = |/3| — \a\, by (11. lip . (12. ip . and Leibniz' formula, for every 
9,a eN'^ we have 

(4.28) |9|9>„,^,,,,„,„...,,^.(x,0l < Cl^°l+l^l+So!5!(l + \x\ + (^^y\<^\ ^ 

for some constant C depending only on 6 and a. 

Since 1 — y9^(x) is supported where |a;| > |/3|^/^, using Leibniz' formula again and 
fliTrjl we get 

X (i + \x\ + \^\r~\'\{xy\-\ 

for some new constant C depending only on 9 and a. We obtain 

(4.29) ||Eo ((1 -y,^(x))p„,;3,.,,o,7i,...,7.(^,^))llB(Q») < C^I+l'l+Sol^ll/^r^. 

Moreover we have 
(4.30) 

1 |«|!|/3 - S\\. f\a\\\P - 5|! V/2 - 5|! , V2 



|a|!V2|/3|!i/2 |«_^|! i^' V |a - 5|!|/3|!y V|«-5|! 

By dUn]), (^M), ^^M, we get 



<1 <1 



(4-31) |^|,i/2|^||l/2 l'^" 

X ||E((1 - ^^(a;))p„,,,,,,„,,,,...,,^,(x,D)(a;'^-^-^°-----^^-9"-^-^°-----^^w))y^ 

p|a| + |/3|-|5|-|(5|-|7o+...+7j|-|7o + ...+7j 
< C fC £:V'^l + l'^l + l'l'0+-+7jl + l70 + ---+7jlr 

- ' ' |/3-5-7o-...-7,|! 

X 1 1 a;^-5-70 - . . . -7, ,9«-5-70 - . . . -7, ^ 1 1 _ 

Since |/3| > |q;|, we have M{a,/3) = |q;|!^/^|/3|!^/^, whereas from fl4.27p we see that 
M{a — 5 — 7o — . . . — 7j, /3 — (5 — 7o — . . . — 7j) = |/3 — 5 — 70 — . . . — 7j|!. Hence we 
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deduce from f l4.3ip that 

X \\E{{1 - M^))p^,^,s,,o,,^,...nA^,mx^-'^^^^^^ 

g|a|+|/3|-|5|-|5|-|7o+... +7,1-170+.. .+7j 



< CsiCsS 



1151 + 151 + 170+. ..+7jl + l70+...+7jL 



M{a - 5 - 7o - . . . - 7j, /3 - 5 - 7o - . . . - 



We now use KW . (I4l9i) . ( lOOl) ( lOHl) . (l432l) to conclude, by the same argu- 
ments as in the case \a\ > that 

(4.33) Yl ]i47^ll^((l-^^(^))[^'^'^1^)llQ=<C':^^^liM- 

\a\ + \l3\<N ^ '' ^ ' 

|/3|>|c,| 

Estimate of ^^ \\E{^p{x)[P,x^d'']u)\\Qs. 

We now start from the formula ( 14 ■41) . For any fixed a,P,S, we choose S < f3 
such that 1/3 — 5| = |a; — 5|, which is possible because here \/3\ > \a\. Writing 
= x^x^~^ in fl4.4p and using (12 .Mp we still get the formula f l4.5p (notice however 
the choice of S is different from the one we made there). We now apply the iterative 
argument detailed at the beginning of the present proof, which led to (14 .Sp , where 
the coefficients Ca,i3,5,'yo,^i,...,-yj satisfy (I4.10p and (14. 9p holds. The multi-indices jj, 
j > 0, are chosen here to satisfy, in addition, = which is possible because 
1/3 — 5| = la — Hence, we can rewrite (14.100 as 

(4.34) < - ' ^2l^o+-+^-^l. 

d!7o!|/3 - d - 7o - . . . - 7j|! 

Now, on the support of we have |x| < 2|/3|^/^; moreover we have \S\ = |/3| — 
|a| Hence it follows from (fLTTj) and ( HTOj) that 

|9|9:(^;3(x)p«,^,,,7o,7„...,7,(x,0)| < Cl^°l+l^l+So!5!|/3|^(l + |x| + |e|)^'"-l^l(x)-H 
for some constant C depending on a, 6. As a consequence, 

(4.35) po (y,^(x)p„,;3A7o,7i,...,7.(^,^))ll6(Q») < C7^l+l^l+So!5!|/3|^. 
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Now we see from Stirling's formula that, for some C > 1, 
(4.36) < C7l/^hH < 

By applying (gSH), fOSj) and (06]) we obtain 

]a\ + \l3\ 

|a|!i7^i/2i 



(4.37) 1^,111/21 fllll/2 l^".^><5.70,7iv,7j I 



5| + |i5| + |7o+...+7j| + |7o + ...+7^ 



|a| + |/3h|5|-|(5|-|7o+...+7j|-|7o + ...+7j 



1/3 - (5 - 7o - . . . -7j|! 

X I U/3-5-70 - . . . -7, 5a-'5-70 - . . . -7, ^ 1 1 



Since > |q;| we have M(q;,/3) = |a;|!^/^|/3|!"'^/^, whereas our choice of 5 and 7^ 
implies that |a; — 5 — 70 — . . . — 7j| = |/3 — 5 — 70 — ... — 7j|. Then we have 
M{a — 6 — % — — 'J j , /3 — 6 — Jo — ... — ■J j) = |/9 — 5 — 70 — . . . — 7j | ! and we can 
rewrite (14.371) as 



< C,(C,£ 



p|a| + |/3|-|5|-|(5|-|7o+...+7j|-|70+---+7j 



||5| + l'5| + l70+---+7jl + l70+...+7iL 



M{a - 5 - 7o - . . . - 7j, /3 - (5 - 7o - . . . - 7j) 

X I U/3-5-70 - . . . -7, ^a-5-70 - . . . -7, ^ 1 1 



It follows from (14. 8 p and (I4.38p . by the same argument^ as in the case |a| > 
that 

(4.39) Yl ]^^7-^P(^/^(^)[^'^'^1^)llQ^<C:e5^-iN- 

\a\ + \fi\<N ^ 
|/3|>|q| 

This estimate, together with (I4.16P and (I4.33p . implies (14. 3p . which concludes the 
proof. □ 



^To be precise, here we do not have longer \S\ < \S\, but rather \S\ = \/3\ — \a\ + \S\, so that 
the number of muhi- indices S which may arise is estimated by 2l'^l~l"l+l^l+''~-'^, and this factor is 
absorbed by the power el''! = el^l-l^l+l^l in (035]). 
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We now turn the attention to the nonhnear term. We first treat the case when 
m > 1. 

Proposition 4.5. Let E e OPr-"(M'^), m > l,h e N,pi,...,p, G N'^, with 
h + max{|pfc|} < m — 1. Let g be a real-analytic function on M.'^ satisfying the 
estimates 

(4.40) |(9"^(x)| < Cl"'+^a!(a;)'^-l"', a; G a G N'^ 

for some C > independent on a. Then for every integer s > d/2 + maxfe{|pfc|} 
there exists a constant Cs > such that, for every e small enough, N > 1 and 
u G iS(M'^), the following estimates hold: 



a\ + \li\ ' 



(4.41) M7^ll^(^'^°(^(^)n^"="))llQ= ^ Cse{S'^_M)'- 

0<|a| + |/3|<Af ^ k=l 

Proof. We first treat the terms with /3 7^ in the sum (14.411) . Let j G {1, . . . , rf} 
such that /3j 7^ 0. By Leibniz' formula, we have 

I I I 

x^d''ig{x)WdP^u)=Xj y V -—^-—d^°g{x)xf'-'^\\d^^^P^u. 

fc=l (5o<a— (5i+...+(5i=Q!— (5o fc=l 

Let now 5o be a multi-index of maximal length among those satisfying |5o| < 
|5o| and 5o < (3 — Cj. Write x^~'^^ = x^°x^~^'~^^ and observe that the symbol 
XjX^'^d^"g{x) belongs to r^"'"^(]R'^), with every seminorm estimated by Al''ol~'"^(5o! 
for some positive constant A independent of 60. Then E o XjX^°d^°g{x) belongs 
to OPr-'"+i+^(M'^). Consequently, it is continuous Q*-^(R"') -> Q"(R"'), with 
M = max{|pfc|}, since —m + 1 + h < —M and its operator norm is bounded by 
yil'5o|+i^pi fQj. j-^g-^ constant A independent of 5o- Then we have 



\E{x^d^{g{x)l[d'''=u))\\Q. 



k=l 



So<a 5i+...+5i=a-5o v ""^ ' ^=1 

We can now write 

I I 



I 



X'" 

k=l k=l 



where 71 + . . .+7; = (3 — ej — 6o and, if \/3\ < \a\, with |7fc| < \Sk\ for 1 < /c < / (which 
is possible because in that case \/3 — Cj — So\ < \a — 6o\; observe that if \So\ < |5o| 
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then /3 — ej — 6o = 0), whereas, if \/3\ > \a\ + 1, with |7fc| > \6k\ for 1 < /c < / (which 
is possible because in that case \So\ = \6o\ and \(3 — ej — 6q\ > \a — 6o\). Moreover, 
if \f3\ < \a\, (then M(a,/3) = |a|!), we have by (12. 3p the following inequality 

la! 1 

(4-42) -jj ■ T < 



a|! 6^1. ..5i\ - M...\6i\V 
whereas, for > |a| + 1, (then M(a,/3) = \a\\'^/^\l3\\^/^), we have 

(4-43) ,„,n/.i.in/. ■ < 



which also follows at once from (I2.3p . Hence by Proposition 12.11 we get 



(4.44) Ji^^\\Eix^d-{9{x)l[d^^u))\\Qs 



5o<a <5i+...+5i=a-<5o k=l ^ 

Let now T G OPr^^''^(]R'^) be any operator which gives an isomorphism Q^~^' — )■ Q'^, 
and write x'^'^d^'-^P'^u = d^'' [x^'^cf'^u) + [x'^'=(9''^ 9^'=]m in the last term of (1131D- We 
get 

\\x^"^d^''^f"^u\\Q.-M < \\x^'^d^^u\\Qs + \\T[x^'d^'^,dP'^]u\\Qs, 

where we used the fact that dP" is bounded Q'(R'^) -> Q"-^^(M'^). 
Using this last estimate we obtain 

f\»\+W\ ' 

' fc=l <5o<« (5o</3-ej- 

(5i + ...+(5;=a-(5o k=l ^ ' l7l<'Ti— 1 

(recall that the 7fc's depend on a, f3,So,Si, . . . ,Si and the choice of Cj). We now 
sum the above expression over |a| + < A^, a 7^ 0. When a and (3 vary but S 
and (5o are fixed, every term in the above sum also appears in the development of 

{ E ^^7^{ll^'^'"ll«^+ E \\T[x^d^d'']u\\Qs]]\ 
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and is repeated at most d times (corresponding to the possibile choices of Cj). 
Hence, taking e sufficiently small, we obtain 

0<|a| + |/3|<JV ^ ^1-^' k=l 

/3#0 



<c''e{ y - — ^wd^u\\Qs+ y mx^d'^.d 

|a| + |/3|<Af-l ^ ' l7l<ni-l 



where we used Proposition 14.41 applied with d'^ and T in place of P and respec- 
tively, and we understand 5'^i[m] = 0. 

We now treat the terms with /3 = in the sum (14.41 1) (recall, M(a,0) = 
Let a 7^ and j E {I, . . . ,d} such that aj ^ 0. By writing 5" = djd^~^^ and using 
Leibniz' formula we have 

k=l 5o<a—ej 5i+...+5i=a~ej—5o k=l 

Observe that Edj o a^o^(x) G OPr-™+i+'^(M'^) is bounded Q'-^^{R'^) Q'(M'^), 
with M = max{|pfc|}, because —m + 1 + h < —M, and its operator norm is 
estimated by A''^°'^^5o! for some positive constant A independent of ^o- Hence 



■\\Edf^{9{x)\[d''^u))\\Q. 



' ' k=l 



8o<a—ej (5i+...+<5;=a— e,— (5o k=l 



Using the inequality 



1 (a-eAl 1 



|«|! S^\...Si\ - M...\Si\V 
Proposition 12.11 and the boundedness of d^'' : Q'^ — )■ Q^^^^ we get 

fe=l (5o<a— (5i+...+(5;=ci— Ej— 5o fc=l 

By the same arguments as above we obtain 

0<la|<Af ' k=l 
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which concludes the proof. □ 
We are now ready to conclude the proof of Theorem 14. 1[ 

End of the proof of Theorem \4.1\ (the case m > 1). From (14. ip we have, for 
a,(3 eN'^,e> 0, 

M+W\ pH+m p\aW\ 

— — -x^cTPu = — — ^x^a"/ + — — 

M(a,/3) M(a,/3) M(a,/3) ^ ^' 

so that 

cri"i+i/3i ^ pH+m „ pH+m „ 



M{a,l3) ^ ^ M{a,[5y ^ M{a,l3) M{a,l3) ^ ^ 

We now apply to both sides the parametrix E of P. With R = EP—I G OPr"^(]R'^) 
we obtain 

M+W\ „ fH+W\ ^ M+W\ 



M(a,/3) M(a,/3) ' ' M(a,/3) 



Taking the Q'^ norms and summing over |a| + < give 

M{a,/3) 



(4.45) S%'[u] < \\u\\qs+ Y1 



0<|Q!| + |/3|<Af 
0<|a| + |^|<iV ^ '^^ 

0<|a| + |^|<iV ^ '^^ 

The second and the third term in the right-hand side of (14.451) can be estimated 
using Propositions 14.31 and 14.41 while the term containing / is obviously dominated 
by S^[f]. For the last term we can apply Proposition 14.51 Hence, we have that, 
for e small enough, 

S'/M < WuWqs + CsS:^V] + Cse{s'^UM + $^(s^l,M)'). 

I 

Iterating the last estimate and possibly shrinking e, we obtain that < oo, 

which implies u e 'HsectiM'^) by Proposition 13.31 
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□ 

4.2. Proof of Theorem I4.lt the case < m < 1. In this case the nonhnearity 
fll.lSp . due to the restriction h+max{\pk\} < max{m— 1, 0} reduces to the foUowing 
form 

(4.46) F[u] = J2Fii^>\ 

I 

the above sum being finite, with I E N, I > 2, and Fi{x) real-analytic functions 
satisfying the following estimates 

(4.47) |9°Fj(x)| < x G a G 

for some C > independent of a. 

We follow the same argument used for the case m > 1, so that we only sketch the 
proof. For technical reasons which will be clear in the sequel, here it is convenient 
to work in the framework of the usual Sobolev spaces, i.e. by defining 
(4.48) 

It is easy to see that the results in Propositions 13.31 14.31 and 14.41 continue to hold 
with and >S'^'^[m] replaced by >S'^^[u] and and with the spaces Q'^ 

replaced by everywhere (operators in OPr'^(]R'^) are bounded on every H'^ by 
Proposition 12.21 with m = n = 0). It remains to estimate the nonlinear term. On 
this point we observe that although this term is more elementary than before, the 
action of the parametrix gives a lower "gain", since < m < 1. Then we have to 
modify slightly our technique. We have the following result. 

Proposition 4.6. Let E G OPr-'^(R"'), < m < 1, and let I e N, I > 2 and g 

be a real- analytic function on satisfying the same estimates as in f l4.47p . Then, 
for every integer s > d/2 there exists a constant > and, for every r > 0, there 
exists Cr > such that, for every e small enough, N > 1 and u & S (M'*) we have 

Proof. We first consider the terms with P ^ 0. We can write 

1 ' 

x'^d''(g(x)u^) = g(x)x'^d''u^ + V -—-^^-—d^'g(x)x^Y\d^^u. 

^-^ On'oi' ■ ■ ■ Oil 

Sq^O 
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Let 60 be a multi- index of maximal length among those satisfying \6o\ < \6o\. 
So < (3. The operators E o g{x) and E o x^°d^°g{x) belong to OPr-™(R'^), and 
the symbol of the second one has each seminorm estimated by Al'^"l'*'^(5o!, for some 
positive constant A independent on 8q. Hence, by the continuity properties on 
weighted Sobolev spaces (Proposition 12.21 with n = 0) we have 

\\E(x^cf(g(x)v})\\„. < CWixyx^x^-'icfiu'^U. 

<So+*l+---+'5fc=« fe=l 

Now, since [3j 7^ for some j e {1, . . . , rf}, we have 

h\<s 

l7l<s 

Now, for every r > there exists > such that 

(4.50) < r|xj| 

Using this inequality and commuting xj with d'^ we get 

5^ ||(x)-'"x,9^(x'^9"(m^))1U2 < rj^ ||a;,9^(x^-^^-9V)|U2 + C,,.||x^-^^-9"(m') 

|7|<s l7|<'5 

<TC,\\x^d-{u^)\\H^ + C:jx^-'W^{u^)\\Hs. 

We notice moreover that for 7' 7^ we have d'^' {{x)'™'Xj) E L°°(M'^), so that 
E E (^,)l|5'''((^)""a;,05^"^'(x^"^^aV)|U2 <C',||X'3-^^■9"M'||H- 
Hence we have obtained that 

{4.51)\\E(x^d''{g(x)u')\\H. < TC.||x''a"('t.')||„. + C;,,||x''-''a°('t.')||„. 

+ E C-l««-^||x-*.n9'.«|U.. 

<S0+il+-..+'5fc=(i ^' k=l 
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Let us estimate the three terms in the right-hand side of (14.511) . To treat the first 
one we observe that 



5i\...Si\ 

(5j + ... + i;=a k = l 

where, as before, we can choose 71, . . . , 7/ G N'^ such that 71 + . . . + 7/ = (3 and 
< (respectively (|7j| > \6j\) if \(3\ < \a\ (respectively if \(3\ > \a\). Then, 
using the same arguments as in the case m > 1, we obtain 

(4.52) rC, -^^J^'d-i^')\\H^<rlCs\M^^^^^ 

\a+0\<N ^ ^ ' 

Similarly, we easily prove that 

(4.53) Y ]^^7^il^'"^^^"(«')lk=<C.^(^^liM)'- 

Concerning the third term in (I4.5ip we can write = nL=i ^^^5 where 71, . . . , 7^ 
satisfy 71 + . . . + 7/ = /3 — 5o and |7j| < \5j\ (respectively > \6j\) if \f3\ < \a\ 
(respectively if \f3\ > \a\). Then, the same arguments as in the case m > 1 yield 
(4.54) 

\c,+f)\<N ^ " So + Si + ...+S^=a ^ ' k=l 

for e > sufficiently small. The estimate of the terms in fl4.49p with /3 = (hence 
a 7^ 0) is very similar but easier, relying on the inequality 

(4.55) (e)-'"lOl<r|Ol+C 

in place of (I4.50p . We omit the details for the sake of brevity. □ 

End of the proof of Theorem 4-i\ (the case < m < 1). Using the same argument 
as in the case m > 1, by the variants of Propositions 13. 3[ 14.31 and with 5'^^[/] 
and S^lf] defined in (14.480 in place of and S^[f], and with the spaces 

replaced by H^, and by Proposition 14.61 we obtain 

S'A^] < \\u\\h^ + + C'.sS'n-A^] + J2 {rC:\\u\\^H}S'/[u] 
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for every N > 1 and e small enough. Now, choosing r < (2 X]/ we 
obtain 

I 

Then we can iterate the last estimate observing that, shrinking r and then e, the 
quantity eCr + r + e can be taken arbitrarily small. This gives < oo and 

therefore u G 'Hsect(^'^)- 

5. Examples and concluding remarks 

5.1. Some remarks on the analyticity estimates. Let us say a few words on 
the estimates 

(5.1) < C'"l+^a!(x)-l"' Va G N"^, a; G R'^ 

which are assumed for the coefficients of the metric in f ll.6p . (cf. also the nonlin- 
earity in dESD, (OD, and f lTT^ ( Kl^ ). 

Locally they are exactly the usual estimates of real-analyticity. To better un- 
derstand the meaning of the decay for |x| — )■ +oo, let us consider the following 
important class of examples. Consider a real-analytic function / in M"^ satisfying, 
in polar coordinates r, w, r > 0, a; G S*^^^, 

f{ru) = h{r~^,u), for r > Tq, u E 

for some tq > 0, where h is an analytic function on [0,r^^) x S'^"^, hence analytic 
up to in the first variable. Let us verify that then / satisfies the estimates f l5.ip . 

Clearly, it is sufficient to check the estimates (15.11) for large Now, by as- 
sumption we have 

oo oo 

/(x) = f{ruj) = E — ¥?fc(w)r-^ = ^k{x), \x\ > Tq, 

k=0 ' k=0 

where ipk are analytic functions on S'^^^, and (pk{x)yipk{'->j)f~^- Observe that the 
functions (pk{x) are real-analytic functions for x 7^ and positively homogeneous 
of degree —k. Moreover, by the very definition of 9?^, for every x, |x| > ro, we have 
the estimates 

(5.2) |9"(^fc(a;)| < C'^l+'^+^a!, 

for some constant C > 0, in some neighborhood of x (this is easily verified in polar 
coordinates and then one uses the analyticity of the change of variables). Hence, 
by compactness, the estimates (15. 2p hold, say, for = 2ro. By homogeneity we 
deduce that 

(5.3) |a"(^fc(x)| <Cl 



X 



2ro 



-K-~\a.\ 



X ^ 0. 
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Hence we obtain 



X 



2roC 



-k—\a\ 

< C(2roC)l"la!|x|-l"l, 



\d^f{x)\<Y,\d''^k{x)\<Ca\Y, 

k=0 k=0 

if |x| > 4roC. This concludes the proof of flS.ip . 

As another remark, we observe that the estimates (15.1 p are in fact equivalent 
to requiring that f{x) extends to a bounded holomorphic function f{x + iy) in a 
sector of the type fll.lSp (see e.g. [HI Proposition 5.1]). This is very useful to check 
the estimates (15. ip in concrete situations, as we will see below. 

5.2. Metric Laplacians. Consider a smooth Riemannian metric gjk{x) in M"'. 
The corresponding Laplace-Beltrami operator has the form 

Cu= Y] r^ ^i {\/g{x)g^^{x)dkU 

-'g[x) ^ > 



j,k=l 



j,k=l \ 1=1 

where g^^ is the inverse matrix of gjk, g = det (gjk), and the Christoffel symbols 
are defined by 

1 

k=l 

Let us assume that the metric is real-analytic and satisfies the estimates 

(5.4) Id'^g.kix)] < gix) > C-\ 

for some C > 0, and every a G N*^, x G W^. Then the matrix g^^ satisfies (II. 5p and 
the estimates in (II. 6p . If in addition we consider V{x) and F[u] as in ( II. 7p . ( II. Sp . 
then the equation 

(5.5) - Cu + V{x)u-\u = F[u], AgC, 

is a special case of (II. 4p . Hence, by Theorem 11.11 every solution u G if*(M'^), 
s > d/2 + 1, of (15. 5p . extends to a holomorphic function u{x + iy) in the sector 
of C"' in (I1.15p . satisfying there the estimates in (I1.16P for some constants C > 0, 
c> 0. 

As a model for the above type of metrics, one may consider the hyperboloid 



5 = {(x, t) G M'^ X M : t = VT+N^}, 

parametrized by x G W^. The Riemannian metric induced on S by the Euclidean 
one then satisfies the estimates ( 15. 4p . as one can easily verify. 
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More generally, we can consider real-analytic scattering metrics in Mf^. They 
play an important role in geometric scattering theory, see [121 EI], P0| Chapter 
6], and have received particular attention in the last years (see e.g. [H] and the 
references therein). Indeed, natural perturbations of the Euclidean metric fall in 
that category. 

A real-analytic metric in M'^ is of scattering type if for any coordinate chart 
V C Ei'^~^ and for some rg > it has the form 

(5.6) h{r^^ ,1]; dr,rdri), for r > Tq, 

where r = |a;|, rj = {rii,ri2, . . . iTjd-i) are real-analytic coordinates on V , and /i is a 
positive definite quadratic form in the last couple of variables, whose coefficients 
are analytic functions on [0,rQ^) x V. Moreover one requires that h{0,rj]dr,dr]) is 
positive-definitejlNotice that this metric approaches the conic metric h{0, rj; dr, rdrj) 
as r — )■ +00, which explains the terminology, sometimes used in the literature, of 
asymptotically conic metric. Notice, by comparison, that the Euclidean metric 
\dx\'^ in polar coordinates reads in fact dr"^ + r^/i', where h' is the usual metric on 

Now, using the remark in Subsection 15. II one sees that, in Euclidan coordinates, 
such a metric satisfies the estimates (15. 4p . In particular, the bound from below in 
(15. 4p is satisfied because h{0,u; dr,rdco) is positive definite, and this last metric in 
Euclidean coordinates has coefficients which are homogeneous functions of degree 
(in fact, each rjj is a real-analytic function of x on ]R_|_ x V, positively homogeneous 
of degree 0). 

5.3. The linear case. This above result for the equation (15. 4p seems interesting 
even in the linear case {F[u] = 0), namely for the eigenfunctions of —C + V{x). 
That equation appears naturally, for example, when looking for standing wave 
solutions (i.e. solutions of the type v{t,x) = e^^^u{x)) of the Schrodinger equation 
idtv — Cv + V{x)v = for scattering metrics (cf. [14j). 

In the linear case we can even assume u G 5'(]R'^). In fact, the existence of 
a parametrix for — £ -|- V{x) (Proposition 12. 3p . shows that such a solution is au- 
tomatically in iS(M'^). Moreover, if V{x) is in addition real- valued, we know e.g. 
from [in] (see also [23] Theorem 4.2.9]) that the operator —C + V{x), regarded 
as a symmetric operator in L'^{R'^,^dx) with domain iS(M'^), is essentially self- 
adjoint and L^(]R'^, ^/gdx) has an orthonormal basis made of eigenfunctionfl Also, 
dimKer (— £ + V{x) — A) < oo, which implies that the width e of the sector in 
(I1.15P can then be chosen uniformly with respect to the solutions. 



Of course, this is equivalent to saying that /i(0, 77; dr, rdrj) is positive definite for every r > 0. 
Since the eigenvalues of the metric are bounded from below and from above, L'^{m.'\dx) = 
L'^(M.'\^dx) as normed space. 



32 



MARCO CAPPIELLO AND FABIO NICOLA 



5.4. Sharpness of Theorem 11.11 We recall that in the case of a differential 
operator with polynomial coefficients the solutions u G iS'(M'^) of the equation 
Pu = extend to entire functions on satisfying estimates (I1.16P in a sector of 
the form fll.l5p . cf. [7J Theorem 1.1]. Very simple examples show that, even in the 
linear case, in Theorem 11.11 we cannot expect an entire extension for the solution 
u. For example, consider, in dimension d = 1, the equation 



where V{x) = + 3 + A solution is given by u{x) = ^^r^e which 

does not extend to an entire function on C. 

The following example shows that, in fact, infinitely many singularities can occur 
along any fixed ray in the complex domain. Let 6 G {—tt/2, vr/2), and consider again 
the equation (15. 7p . with 



By applying the remark at the end of Subsection 15.11 to the function tanh(e~*^a;) 
it is immediate to check that V{x) satisfies the estimates in (11.71) . On the other 
hand, the function 



is a solution of (15. 7p and extends to a meromorphic function in the complex plane 
with poles at ;z = e^(^+^/2)(2A; + l)7r, k E Z. 

This shows that in Theorem II. H even in the linear case, the form of the domain 
of holomorphic extension as a sector is sharp in general. The following example 
shows a similar phenomenon in the nonlinear case, even for the standard harmonic 
oscillator. 

Consider the following nonlinear perturbation of the harmonic oscillator, in di- 
mension d = 1, ai the first eigenvalue A = 1: 



(5.7) 



u" + V{x)u = 0, 



V{x) =x^ -I- 2e-'^'^ + 4e-*^xtanh(e-*^x) + Qe-^''^xia.nh^{e-'^x). 



u{x) = cosh-'(e-'^x)e-"'/2 



(5.8) 




It was shown in [7] that the solution of (15.80 is given by 



(5.9) 




with A = M| 
by 



,1-fc 

'0 



, where we used the complementary error function defined 



Erfc(t) 
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Here and in the following, roots are defined to be positive for positive numbers, 
with continuous extension to the complex domain, i.e., we take principal branches. 

Suppose now A > 0, that is < Mo < (^^) • In this case, smce 

< A < A + V2A;-2Erfc(^^^^a;^ , 

the solution u{x) in (15. 9p is well defined analytic in M and 

< u{x) < Ai-^e 2 . 

Similar estimates are valid for u'{x), u"{x). Hence we have u G if^(M), and The- 
orem 11.11 applies, implying the desired holomorphic extension u{z) to a sector. 
However, as observed in [7], u{z) is not entire, but has a singularity dX zq E'C when 

(5.10) A + ^2k - 2Erfc(^y^^Zo^ = 0, 

where Erfc(z) is the entire extension of Erfc(x). Such singularities in fact occur, 
because the great Picard theorem in the complex domain grants the existence of 
infinitely many solutions of (I5.10p for all A G C, but for a possible exceptional 
value, see [29] . 

Indeed, we now prove the following more precise result. 

Proposition 5.1. For every A > 0, hut for a possible exceptional value, and every 
e > 0, u{z) has a sequence of singularities which tends to infinity in the sector 
7r/4 < arg^; < 7r/4 + e or in 3'7r/4 — e < aigz < 3tt/A. 

Proof. Using the great Picard theorem as above and the refiection properties 



Erfc(z) = Erfc (2), Erfc(-z) = y/^ - Erfc(^) = - Erfc(z), 
which can be verified directly from the definition, it is sufficient to prove that 

(5.11) Erfc(2;) -H- as 2; — )■ 00 in the sector |arg2;| < 7r/4, 
and that, for every e > 0, 

(5.12) |Erfc(2;)| — +00 as 2; — )■ 00 in the sector tt/A + e < a.Tgz < 7t/2. 
Now, (15. lip follows at once from the expansion 



,2 



Erfc(z) = ^ — {1 + R{z)) with \R{z)\ < ^ 



2z ' ^ " ' ' " - V2I 

valid when |arg2;| < 7r/4; see e.g. [HI pages 18-20]. 
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The property f l5.12p can be verified directly as follows. Observe that, for z = 
X + iy, X > 0, y > 0, we can write 

Erfc(2) = - f e-'^'du, 

J 'y 

where the path 7 is given by the hyperbola through z = x+iy with parametrization 
u = 7(t) = f+it,te (0, y]. Then 

\ETic{x + ty)\ = + 



Let < ;U < 1 be a number to be chosen later. We have 

|Erfc(x + «y)| > r e-"^^'' dt > (1 - /i)i/e-'^"'"'+^'^'. 

Jliy 

Now, if z belongs in addition to the sector in fl5.12p . we have < x < ey, for some 
e < 1. We obtain then 



|Erfc(x + iy)| > {1 ~ fx)ye^>' . 

If we choose // > y/i, we get |Erfc(x + iy)\ — > +c>o as ?/ — )■ +00, which gives the 
desired conclusion when x = Re z > 0. The case when x = is immediate, because 

ry 

Erfc(i?/) = - / e*^ dt 



2 

Property f l5.12p is then proved. □ 
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